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Abstract. The thermal efficiency of the kinesin cycle at stalling is presently a matter of some debate, with 
published predictions ranging from (A. W. C. Lau, D. Lacoste and K. Mallick, Phys. Rev. Lett. 99, 
158102 (2007); D. Lacoste, A. W. C. Lau and K. Mallick, Phys. Rev. E78, 011915 (2008)) to 100% (G. 
Oster and H. Wang, in Molecular Motors, edited by M. Schliwa (Wiley- VCH Verlag GmbH, Weinheim 
(2003), p. 207). In this note we attemp to clarify the issues involved. We also find an upper bound on the 
kinesin efficieny by constructing an ideal kmesin cycle to which the real cycle may be compared. The ideal 
cycle has a thermal efficiency of less than one, and the real one is less efficient than the ideal one always, 
in compliance with Carnot's theorem. 

PACS. 87.16.Nn Motor proteins - 05.40.-a Fluctuation phenomena, random processes, noise, and Brow- 
nian motion 

The application of physics thinking to life is one of The above said is particularly true of the application of 
the great breakthroughs of modern science [IJ. The big thermodynamics to biophysical problems 5J. If a physical 
questions (what is the fundamental difference between liv- description of a living system is possible at all, it must in- 
ing and nonliving matter? is there life anywhere else in elude thermodynamics at some point. But when we focus 
the Universe? [5]) are conspicuously evading us, but the on phenomena at the molecular level, it is the thermo- 
progress in the analysis of concrete phenomena is aston- dynamics of small systems we are talking about, not the 
ishing [Sill] . Since this generally involves applying physical familiar macroscopic thermodynamics of heat engines and 
theory in contexts quite removed from the original ones, stars [B]. 
we physicists are learning quite a lot in the process. 

The analysis of molecular motors is a case in point. 
Analyzing engines is what thermodynamics is good for; it 
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was invented for just this purpose. But when we get to a 
particular molecular motor the possibility (let alone the 
usefulness) of a thermodynamic analysis seems somewhat 
less than obvious, and one correspondingly finds a variety 
of points of view and statements in the literature, even 
when addressing the same systems. 

To be concrete, let us take the example of kinesin re- 
garded as a molecular motor. Kinesin converts chemical 
energy into useful work; it should be possible to define 
the thermal efficiency as the ratio of output work to input 
energy. One expects this should be less than one by a fi- 
nite amount, lest some kind of perpetual mobile could be 
built on kinesin. Indeed one can build a model of kinesin 
which predicts that the efficiency drops to at stalling 
[7ll8j . However a careful analysis reveals that the opposite 
conclusion of 100% efficiency also has some grounds [5]. 

Let us point out that the problem is not that there are 
several nonequivalent measures of efficiency going around 
[TD]. We are talking about thermal efficiency [TT1|12| . and 
that is usually clearly distinguished from alternatives such 
as Stokes efficiency ^T3\ . 

To understand how kinesin could reach the 100% effi- 
ciency mark it is instructive to compare the kinesin cycle 
to any well known thermal engine, such as a gas engine 
going through an Otto cycle. The comparison is not arbi- 
trary, since both kinesin and the gas engine extract their 
power from a chemical reaction, and are in contact with 
a single environment [14 . However, there is a potential 
difference between both machines. In the car engine, the 
fuel delivers energy to the working fluid under the form of 
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heat. Therefore the entropy of the working fluid increases 
during combustion. Because of Clausius' inequality, the 
rest of the cycle cannot be completed adiabatically. To el- 
liminate the surplus entropy, heat must be lost to the en- 
vironment, and this detracts from the engine's efficiency. 
Note that this is true even under ideal conditions. After 
every form of friction is elliminatcd, efficiency remains less 
than one by a finite amount. Moreover, the ideal Otto is 
less efficient than the ideal Carnot cycle, if we set the hot 
and cold sources of the latter at the extreme temperatures 
of the former. 

In the case of kinesin, however, it is not obvious that 
heat is involved at any step in the process. Kinesin extracts 
its power from the free energy liberated by the hydrolysis 
of one ATP molecule. One can imagine that this free en- 
ergy is delivered under the form of work and it is stored 
at intermediate steps as potential energy, for example, as 
elastic or electrostatic energy in the kinesin macromolecule 
|15j . Work exchange is free from the restrictions imposed 
by the Second Law on heat exchange and heat to work 
conversion, and so a 100% efficiency is not excluded. 

Therefore the issue hinges on whether at any step within 
the kinesin cycle we may identify a conversion of work into 
heat, or at least an entropy increase. If this is the case, 
then the kinesin cycle cannot be completely adiabatic, and 
the efficiency must be strictly less than 100%. Moreover, 
this entropy increase must obtain even after known forms 
of friction, or of entropy production because of processes 
evolving at a finite rate, are elliminated. In other words, 
we ought to be able to identify an ideal kinesin cycle, and 
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show that under no conditions the ideal kinesin cycle gets The rest of the paper is organized as follows. In the 
to 100% efficiency. next section we provide some basic background material 

on the kinesin cycle. Then in the following section we 

In a recent contribution, Bier [16j has presented a sce- 

present the Crooks fluctuation theorem and we argue that 

nario along these lines. The kinesin cycle is an asymmetric 

this theorem can be applied to the kinesin cycle. In Sec- 
hand on hand walk along a microtubule. At some point in 

tion IV we use it to estimate the free energy change and 

the kinesin cycle, both kinesin heads are attached to the 

the efficiency of the process. We conclude with some very 

microtubule. Then one of the head detaches. According 

brief final remarks. 

to [16] , while the head of kinesin is free it has many more 
available states than when docked, and so it also has a 

1 A brief description of the kinesin cycle 

larger entropy. Therefore docking must be accompanied 

by an entropy decrease, and a corresponding heat loss Q Let us begin by stating the essential facts about the ki- 
to the environment. The kinesin efficiency must therefore nesin cycle pOll^l^^ll^ . 

be less than 100%. Kinesin is a polymer which can attach to a load on one 

end and has two heads on the other end. By attaching the 

Our goal is to further advance this argument by es- 

heads to a microtubule, kinesin can perform an asymmet- 

timating the free energy change AJ- associated with a 

ric hand on hand one-dimensional walk, effectively drag- 
forward step of kinesin along the microtubule. Because ki- 

ging the load against an external force F and the viscous 

nesin is in contact with an environment at constant tem- 

drag from the environment, thus doing useful work. Each 

perature, the free energy change poses an upper bound 

step in this walk corresponds to a cycle. The strides are 

on the work W which may be extracted from kinesin, 

S — 8 nm long, corresponding to the distance between 

W < {—AT), with equality obtaining under ideal con- 

the attachment sites on the microtubule. During one step 

ditions. 

one head remains firmly docked, while the other head tra- 
The problem is how to extract AT from actual exper- verses double that distance. We give in figures [1] to |4| a 
iments, since these are done under non ideal conditions, simplified account of the kinesin cycle. 
To overcome this difficulty we note that, if we regard a One crucial input for the analysis below is the amount 
backward step as the time-reversal of a forward step, then of energy liberated by the hydrolysis of one ATP molecule. 
AT can be related to the ratio of the probabilities for for- We adopt the estimate provided by Bier [16], namely 25 ksT. 
ward and backward steps through the Crooks fiuctuation Here ks = 1.4 10^^"^ J K^^ is the Boltzmann constant 
theorem. Since these probabilities can be measured [171 and T = 300 K is the environment temperature. Let us 
[T51ll9j . this solves our problem. observe that a lower value would change the numerical 
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Fig. 1. At tire beginning of the kinesin cycle, both heads are 
attached to the microtubule. The rear head is associated to 
a molecule of adenosine diphosphate (ADP). The front head 
captures a molecule of adenosine triphosphate (ATP) from the 
environment. 




Fig. 2. In the second step of the cycle, the rear head separates 
from the microtubule and the whole molecule swings around 
the front head, bringing the other head close to the next dock- 
ing site. 




Fig. 3. In the third step of the cycle, the now front head 
releases its ADP molecule into the environment, while in the 
other head ATP is decomposed into ADP and a phosphate 
(Pi), which is released. This process, namely the hydrolysis of 
ATP, is the source of the energy for the whole cycle. 




Fig. 4. Finally, the front head binds to the microtubule. The 
cell disposes of the released ADP and Pi and provides a new 
ATP molecule. This is analogous to the exhaust stroke at the 
end of the Otto cycle. 

prediction for the efficiency, but not the essentials of the 
arguments below; only a value as low as half of this would 
validate the prediction of 100% efficiency at stalling un- 
der ideal conditions. We assume that exactly one ATP 
molecule is decomposed in each cycle, so the number r 
of steps per unit time is equal to the reaction rate for 
ATP hydrolysis. This reaction rate is regulated by the cell 
through the concentrations of the relevant substances. The 
dependence of r with the load is weaker than other rele- 
vant factors [12] and we shall disregard it. 

Upon release, kinesin undergoes a configuration change 
which pushes the free head to the neighborhood of next 
unoccupied docking site. Kinesin can exert a constant 
force, mostly independent of the opposing external force. 
Otherwise, the free head is caught up in the thermal Brow- 
nian motion generated by the environment. Therefore the 
end of the cycle is undetermined: the free head can dock 
in the intended site, in which case the step is successfully 
completed, or else it may be dragged back to its initial po- 
sition. Since the reaction rate is fairly constant, the ratio 
of success to failure is the main determinant of the overall 
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progress rate. The probability of a successful forward step 
over that of failure has been measured to be [l9lfT6] 

where Fgt — 7 pN is the stalling force. At higher loads, 
kinesin drifts backwards. 

It follows from eq. ([T]) that the maximum work kinesin 
can do against the external force is FstS — 13.3 /csT, 
which is about half of the available energy. This puts an 
absolute limit on thermal efficiency at about 50% at any 
operating condition, given our previous estimate for the 
available energy. 

In the real world, there are many sources of dissipation 
that help to explain the missing efficiency, such as viscous 
drag or friction related to the conformational changes in 
the polymer. However, in devising an ideal cycle we are 
free to imagine that those sources of dissipation may be 
controlled. For example, we may dream of a reversible ki- 
nesin cycle where the free head drifts so slowly that it is 
always in thermal equilibrium with the environment. In 
this case, the net viscous drag would be negligible. 

2 The Crooks fluctuation theorem 

When the kinesin cycle is completed, the only changes 
are that the molecule has been displaced along the micro- 
tubule and work has been done against an external force, 
if any. Since the whole process is isothermal, the work W 
extracted is bounded by the free energy drop when kinesin 
moves one step forward. The efficiency of the ideal kinesin 
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cycle, when equality obtains, is therefore equal to the ratio 
of the free energy drop to the input energy. 

To obtain the free energy change associated to one step 
from observable data, we shall apply the Crooks fluctua- 
tion theorem f241l251l26l[?7ll281l29j to kinesin. Let us first 
make explicit our conventions. We call Q the heat taken 
in by the system in a nonequilibrium process, and W the 
work given out by the system. The First Law is written 
a.s AU = Q — W, where U is the internal energy, and the 
Clausius inequality is Q < TAS, where T is the temper- 
ature and S is the entropy. We define the free energy as 
T = U—TS, so the Clausius inequality implies W < —AT 
in any isothermal process. For present purposes, we do not 
need to discriminate between internal energy and enthalpy 
H = U + PV, and therefore we do not make a distinction 
between the (Helmholtz) free energy T and the Gibbs free 
energy g = H -TS. 

The Crooks fluctuation theorem concerns a system in 
interaction with an environment at temperature T. The 
system depends on some external parameter A. Initially 
the system is at state A. The parameter A is made to 
depend on time following a trajectory A/. We call Pf the 
probability that the system ends up in state B, giving out 
work W. If now the evolution of A is reversed, we call P^ 
the probability that the system, now starting from B, ends 
up in state A giving out work — W. The Crooks fluctuation 
theorem states that 
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There are several proofs of the Crooks fluctuation theo- 
rem in the hterature, each based on different assumptions 
on the microscopic dynamics of the system in question. 
Since we lack a microscopic model of kinesin, none of those 
proofs can settle the question of whether the Crooks fluc- 
tuation theorem can be applied to the kinesin cycle. 

A relevant observation is that the kinesin molecule is 
a very complex system but only a few of the variables 
defining its state (such as the position of the heads and 
whether they are associated to ADP and ATP molecules) 
are relevant to our discussion. While formal proofs of the 
Crooks fluctuation theorem usually assume that the initial 
and final states are equilibrium states, from the physical 
point of view we may expect that the theorem holds if the 
characteristic times of the process are long enough that 
the unobservable part of the theory (the particular micro- 
scopic states of system and environment) is capable to ad- 
just to the instantaneous value of the relevant macroscopic 
variables. This only puts a mild bound on the character- 
istic times associated to the macroscopic evolution of the 
system. To make this bound precise we need an estimate 
for the relaxation times of the microscopic model, which 
we do not have at present. 

This expectation is consistent with related work on 
fluctuation theorems applied to systems in nonequilibrium 
steady states [501151] and also on periodically driven sys- 
tems [5 ^ 155 ] . 

The above argument only makes the wide applicability 
of Crooks fluctuation theorem plausible, but it is not con- 
clusive by itself. In last analysis, we must consider whether 
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the Crooks fluctuation theorem has been applied to situ- 
ations comparable to the kinesin process. 

Probably the most detailed match of the Crooks fluc- 
tuation theorem against experimental data is the one re- 
ported in 28J . In this opportunity, the Crooks fluctuation 
theorem was used to estimate the free energy difference 
associated to the unfolding of a RNA molecule. In the 
experiment, a single molecule was repeatedly folded and 
unfolded. We find the periodic folding of a single molecule 
is analogous to the cycles of the kinesin displacement, and 
therefore the success of the application of the Crooks fluc- 
tuation theorem in one case argues for the feasibility of its 
application in the other. 

3 Kinesin and the Crooks fluctuation theorem 

To apply the Crooks fluctuation theorem to kinesin, we 
identify the macroscopic initial and final states A and B 
as the initial and final states in the kinesin cycle, and the 
external parameter A as any variable that measures the 
progress of the molecule from one pair of docking sites to 
the next, such as the expectation value of the position of 
the center of mass of the molecule along the microtubule. 
We also regard a backward step as the time reversal of a 
forward step, both under the same external force F (so 
that work is reversed). Noting that W = F6 for a forward 
step, the Crooks fluctuation theorem [2] with the probabil- 
ity ratio given in [T] implies 

AJ^=^[Fst+F] (3) 
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On the other hand, ~AJ- is the maximum work which may 

be extracted from the system at constant temperature. {W) — F (5) (6) 

We have given above an estimate of the energy input as , ic\ ■ ,i i- i 

where {o} is the average displacement 

25 ksT. This is close to 25Fst = 26.6 kBT[16\. Adopting 

this approximate energy input, the efficiency of the ideal _ tJtanhi ^ [F F]\ (7) 



kinesin cycle is 



1 

V ideal ^ 



and the actual efficiency is 



(4) F (S) 

- 2^ 



Under ideal circumstances, all the energy 25Fst given to 
the polymer at the beginning of the cycle goes either into 
the reversible work — AT or into the heat dissipated dur- 
ing docking. If docking is isothermal, then we have 



Sfree " Sdock = = [SF^t - F] (5) 



Observe that rj^eai vanishes at both gliding and stalling, 
in agreement with [7,8 . In general, we have rjreai < Videai 
always, also in agreement with Carnot's theorem: no ac- 
tual machine can perform better than a reversible machine 
with the same energy input. 

4 Final Remarks 



where Sfree and Sdock are the entropies of the free head 

and the docked states, respectively. The issue of the efhciency of molecular motors, with ki- 

We see that the efficiency of the ideal kinesin cycle nesin as a particular case, is a matter of debate in the 

ranges from 25% when gliding to 50% at stalling (of course, literature. If one regards the kinesin cycle as a series of 

the second value will change if the assumed input energy conversions of one form of potential energy into another, 

is modified). No conflict with classical thermodynamics always through the action of work, then the efficiency can 

is apparent. The efficiency of the real cycle is less than be very high. In particular, at stalling, when no system- 

this. There are several reasons for this reduced efficiency, atic viscous drag opposes motion, the efficiency could be 

of which two are outstanding. First, not all the available as high as 100%. 

work is actually extracted: some of it is left to be dissi- On the other hand, if one adopts the phenomenological 

pated as heat, for example by opposing the viscous drag definition of efficiency eq. [51 then one can say that the 

or as excess kinetic energy to be absorbed by the docking efficiency is zero at stalling by definition, 
site. Second, the cycle may fail, with the polymer drift- In this note, we have attempted to clarify the issues 

ing backwards rather than forward. Therefore the actual involved by noting that for a motor like kinesin, operating 

average work is at constant temperature, the efficiency is limited by the 
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free energy change associated to each step forward along 5 Acknowledgments 
the microtubule. 



The most important contribution of this note is the 
suggestion that this free energy change can be extracted 
from actual experimental data on the ratio of probabilities 
for a forward and a backward step, through the application 
of Crooks fluctuation theorem. Our actual analysis yields 
an efficiency of 50% at stalling. Therefore the prediction 
of 100% efficiency ought to be discarded. On the other 
hand, our analysis concerns only an ideal situation, and 
therefore does not conflict with a prediction of 0% under 
physical circunstances. 

The central issue is therefore if the application of Crooks 
fluctuation theorem is warranted. It is possible to verify 
the Crooks fluctuation theorem for a given microscopic 
dynamical model, but not to actually prove its validity 
in an experimental situation. Very much the same could 
be said about the Second Law itself. In this case, support 
for the application of the Crooks fluctuation theorem to 
kinesin may be found by noting the extreme generality of 
this statement, and that it has already been applied in 
similar situations. 

The ultimate motivation of our work is that the ther- 
modynamic analysis provides guidance for further micro- 
scopic studies of kinesin as a dynamical system. From this 
point of view, we expect the final test of the hypothesis 
advanced in this note will be its heuristic value. 
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